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Abstract. We present results for a numerical study of the ghost propagator in Coulomb gauge whereby lattice results for the 
spatial gluon propagator are used as input to solving the ghost Dyson-Schwinger equation. We show that in order to solve 
completely, the ghost equation must be supplemented by a boundary condition (the value of the inverse ghost propagator 
dressing function at zero momentum) which determines if the solution is critical (zero value for the boundary condition) or 
subcritical (finite value). The various solutions exhibit a characteristic behavior where all curves follow the same (critical) 
solution when going from high to low momenta until 'forced' to freeze out in the infrared to the value of the boundary 
condition. The boundary condition can be interpreted in terms of the Gribov gauge-fixing ambiguity; we also demonstrate that 
this is not connected to the renormalization. Further, the connection to the temporal gluon propagator and the infrared slavery 
picture of confinement is discussed. 
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The ghost propagator and its Dyson-Schwinger equa- 
tion [DSE] are of central interest to studies of nonper- 
turbative Yang-Mills theory. In Coulomb gauge, there is 
known to be two types of infrared behavior: critical (with 
an infrared divergent dressing function) and subcritical 
(infrared finite dressing) [1], similar to the situation in 
Landau gauge [2]. We demonstrate here how the exis- 
tence of these different solutions is intrinsic to the stan- 
dard ghost DSE in Coulomb gauge and how a particu- 
lar pattern emerges. This talk is based on the work of 
Ref. [3]. 

To solve the ghost DSE, we require two pieces of in- 
put: an Ansatz for the ghost-gluon vertex and the spatial 
gluon propagator. In Coulomb gauge, the ghost-gluon 
vertex renormalization coefficient is nonperturbatively fi- 
nite and since the vertex reduces to its tree-level form 
as one ghost leg momentum vanishes [4], we can take 
Za = 1 • Our vertex Ansatz is then to simply take the tree- 
level result (see Ref. [3] for details of our conventions) 
as is standard: 
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The spatial gluon propagator is decomposed as follows: 
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where tn is the transverse spatial projector and the dress- 
ing function reflects the noncovariance of the gauge. As 
it turns out, the important quantity is not Daa itself, but 
rather the dressing for the equaltime propagator 
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For definiteness we use the lattice results of [5], compat- 
ible with a Gribov formula [6], as input: 
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Many other studies have reported similar results (see 
Ref. [3] and references therein for a discussion), moti- 
vating the usage here. The above form is already renor- 
malized and it was found that the Gribov scale m is inde- 
pendent of the renormalization scale jj,. 
The unrenormalized ghost DSE reads [7]: 
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The ghost propagator (W a ) and the ghost proper two- 
point function (r ?c ) are decomposed as follows: 

W*(k) = -8" b ^-D(k 2 ), lf c {k) = 8 ab ik 2 T{k 2 ). (6) 
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The dressing functions obey D{k 2 )T{k 2 ) = 1. Notice that 
as a general result of the Slavnov-Taylor identities [8], 
the ghost dressing functions are strictly independent of 
the energy. With the vertex Ansatz, Eq. (1), one thus 
recognizes the temporal integral, Eq. (3), which can be 
replaced by the Gribov formula, Eq. (4). Introducing 
the ghost propagator renormalization coefficient Z c and 
writing in terms of the ghost dressing functions, the 



renormalized ghost DSE now reads 
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where all quantities are implicitly renormalized at some 
scale jj. and we have N c colors. 

At one-loop perturbatively (obtained by setting m = 
0), using dimensional regularization and the renormal- 
ization scale ji, one obtains the standard result [7] 



r(je) = l + -AA/ r c ln(x/M) 
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where x = k? and X — g 2 / (An) 2 . One can see that for 
small x, T(x) could become negative. Gribov's original 
derivation of the formula, Eq. (4), was based on expressly 
prohibiting such a situation [6]. We shall return to this 
topic later. Naively resumming the result, Eq. (8), gives 
the ultraviolet [UV] asymptotic behavior for the ghost: 
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where y g is (naively) the leading order expression for the 
ghost anomalous dimension. 

One can also investigate the infrared [IR] behavior 
of the ghost DSE. Using a UV-cutoff regularization and 
subtracting at ji, the DSE has the form (/ is a function 
arising from the angular integration, see Ref. [3] for the 
explicit form) 

T(x) = T(p)-XN c J o [I(x,y;m)-I(ii,y;m)} . 

(10) 

The low x and large A scales separate because Za = 1. 
Assuming a powerlaw behavior F(x) ~ (x/ji) a as x — >• 0, 
one obtains the relationship 
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This type of relationship is quite familiar in Landau 
gauge studies (e.g., [9]). Counting powers of x, one sees 
that if there is an IR powerlaw, it will have the exponent 
a = 1/2. However, counting powers of fi instead, one 
would have a = which indicates an IR finite dressing 
function (and is intimately connected to the presence of 
the additional scale m within the function /). 

Equation (10) is not particularly suitable for numerical 
analysis because (as noted perturbatively) the subtraction 
at jx might lead to a negative result and numerical insta- 
bilities. Taking inspiration from Gribov [6] and demand- 
ing that T{x) be positive for all x > 0, one can subtract the 
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FIGURE 1. Results for the ghost propagator dressing func- 
tion D(x). See text for details. 



equation at x = 0. Further scaling all dimensionful quan- 
tities with m (e.g., x — > xm 2 ), the ghost DSE can then be 
written 
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To evaluate this equation, there are two input parameters: 
XN C (X = a s /(4n) with a s = 0.1187 and N c = 3) and 
r(0) > 0. It is clear that T(0) is not given by the DSE, 
but is a boundary condition for the solution. Results for 
D(x) = r(x) _1 are shown in Fig. 1 for various values of 
r(0). In the UV, the solutions with T(0) < 0.3 lie on top 
of each other and agree asymptotically with the naively 
resummed perturbative result, Eq. (9). Going from high 
to low momenta, there is a characteristic behavior where 
the solutions for successively lower values of T(0) fol- 
low the same curve further into the IR until 'freezing' to 
the value r(0) required by the boundary condition. As 
r(0) — > 0, one sees the emergence of the critical solu- 
tion with the powerlaw exponent a = 1/2 as predicted. 
Lattice results for the ghost propagator seem to rule out 
r(0)>0.2[3]. 

It is interesting to also calculate the renormalization 
coefficient for the ghost, Z c . In terms of the scaled vari- 
ables, this is given by 
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Numerically, one can verify that Z c is independent of the 
momentum x, confirming the multiplicative renormaliz- 
ability of the ghost DSE here. Seeing the explicit depen- 
dence on r(0) in the expression above, one can plot Z c 
as a function of r(0) (with fixed A). This is shown in 
Fig. 2. Surprisingly, for r(0) < 0.3, Z c shows almost no 
variation and indicating that the boundary condition is 
not connected to the renormalization. 
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FIGURE 2. Results for the ghost renormalization coefficient 
as a function of T(0). See text for details. 



After eliminating the possibility that T(0) is connected 
to the renormalization, the characteristic pattern of the 
results exhibited in Fig. 1 can be naturally interpreted 
in terms of Gribov's gauge fixing ambiguity [6] (indeed, 
this ambiguity was the origin of the condition T(0) > 0). 
Central to the gauge fixing is the Faddeev-Popov [FP] 
operator, 
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and the Gribov problem arises when zero modes of the 
FP operator occur within the functional integral. (Resolv- 
ing the spatially independent zero modes that are par- 
ticular to Coulomb gauge leads to a conserved and van- 
ishing total color charge [10].) Perturb atively, when one 
expands around A = 0, there are no zero modes but non- 
perturbatively in the IR, A can become large and zero 
modes appear. With zero modes present, the gauge is not 
completely fixed and the functional integral contains an 
integration over the gauge group which when taking the 
expectation value of a gauge dependent quantity would 
lead to a suppression. In this case, the suppression would 
presumably be a freezing at nonperturbative scales. The 
DSE (at least in standard form) on the other hand does 
not change its form in the presence of the Gribov am- 
biguity - it merely relates the Green's functions to one 
another with their exact definition left implicit. Thus, the 
common perturbative curve and the freezing of the so- 
lutions in the IR for finite T(0) would appear to signal 
the presence of the Gribov ambiguity within the implicit 
definition of the functional integral for the ghost propa- 
gator. Conversely, the T(0) = (critical) solution would 
appear to correspond to a more proper definition of the 
ghost. Such a situation can be explicitly shown in 1 + 1- 
dimensional Coulomb gauge [11]. 

The critical ghost solution can be tentatively linked to 
quark confinement [3]. The Slavnov-Taylor identities [8] 
tell us that the dressing for the temporal gluon propagator 
has the form 



where is the longitudinal part of the spatial gluon 
polarization. The renormalized DSE for Faa would 
schematically be 
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where each self-energy loop E, is multiplied by a ver- 
tex renormalization coefficient. After subtracting at the 
scale fx, this would be similar to the form of the ghost 
DSE, Eq. (10), but with nontrivial Z, factors - only the 
ghost-gluon vertex (for which Za = 1) has a trivial co- 
efficient but since the corresponding loop involves en- 
ergy independent ghosts, this loop must cancel (the in- 
famous energy divergence problem of Coulomb gauge). 
Thus, the separation of the external momentum and UV- 
cutoff scales necessary for the IR powerlaw analysis of 
Eq. (10) is no longer applicable and one is left with the 
conclusion that Faa must surely be a constant in the IR. 
See Ref. [3] for other arguments. So with a critical ghost 
solution, Dqq <~ 1/k 2 as k 2 — > 0. This is exactly the be- 
havior one would expect in order to have a linear rising 
potential between quarks in Coulomb gauge [12]. 
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